Observation of squeezed states with strong photon number oscillations 
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Squeezed states of light constitute an important nonclassical resource in the field of high-precision measure- 
ments, e.g. gravitational wave detection, as well as in the field of quantum information, e.g. for teleportation, 
quantum cryptography, and distribution of entanglement in quantum computation networks. Strong squeezing 
in combination with high purity, high bandwidth and high spatial mode quality is desirable in order to achieve 

significantly improved performances contrasting any classical protocols. Here we report on the observation of 
the strongest squeezing to date of 1 1 .5 dB, together with unprecedented high state purity corresponding to a vac- 
uum contribution of less than 5%, and a squeezing bandwidth of about 170 MHz. The analysis of our squeezed 
states reveals a significant production of higher-order pairs of quantum-correlated photons, and the existence of 
strong photon number oscillations. 



I. INTRODUCTION 

Squeezed states as well as number states (Fock states) are 
so-called nonclassical states. They allow the measurement, 
the communication and the processing of information in a way 
not possible with coherent states that are governed by vacuum 
fluctuations. Squeezed states of light and correlated photon 
pairs have been used in order to realize interferometric mea- 
surements with sensitivities beyond the photon counting noise 
[1, 2, 3, 4, 5, 6], to demonstrate the Einstein-Podolski-Rosen 
paradox [7, 8, 9, 10], as a resource for quantum teleportation 
[11, 12, 13] and for the generation of Schrodinger kitten states 
for quantum information networks [14, 15]. In view of appli- 
cations in long distance quantum communication, purification 
and distillation of entangled photon number states and entan- 
gled two-mode squeezed states were experimentally demon- 
strated [16, 17, 18J. 

Fock states are characterized by photon counting detectors, 
whereas squeezed states are traditionally characterized in 
the phase space of position and momentum-like operators 
[19, 20]. Appropriate operators are the non-commuting am- 
plitude and phase quadratures of the field, Xx and X2, respec- 
tively. Their variances obey a Heisenberg uncertainty relation, 
A^Xi ■ A^X2 > 1/16, with the vacuum noise variance nor- 
malized to 1/4. A squeezed state is realized if the quantum 
noise of one quadrature is squeezed below the vacuum noise 
level, at the expense of increased noise of the other quadrature. 
For an overview we refer to Ref. [21]. The relation between 
the quadrature operators and the photon number operator n is 
revealed by the Hamilton operator of the quantized harmonic 
oscillator: 

H = hui{h+ 1/2) = nw{Xl + Xl) , (1) 

where h is the reduced Planck constant and lo the light's angu- 
lar frequency. One can think of the two expressions in equa- 
tion (1) corresponding to the particle picture and the wave pic- 
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ture of quantum physics respectively. These can be used to 
calculate the mean photon number of squeezed states. Note, 
that this number is always greater than zero. Consider now 
a squeezed state without displacement in phase-space. Such 
a state is called a squeezed vacuum state, and one can write 
{Xf) = A^Xj. If such a state is pure, i.e. has minimum un- 
certainty, only even photon numbers can be observed [22, 23]. 
Furthermore, if the squeezing effect is strong a so-called odd- 
even photon number oscillation is realized, that includes not 
just the existence of photon pairs but also the existence of 
higher even-numbers of quantum-correlated photons. 
This article presents the generation and characterization of 
squeezed vacuum states with strong odd-even photon nimi- 
ber oscillations revealing a significant contribution of corre- 
lated photon numbers up to 10 and higher. This result is made 
possible because our squeezed vacuum states show unprece- 
dented quaUty, i.e. are strongly squeezed, are of high purity, 
and are in a well-defined spatial Gaussian fundamental mode. 
These properties are of great importance for apphcations in 
high-precision measurements as well as in quantum informa- 
tion. To the best of our knowledge, our squeezed light source 
also shows the broadest squeezing bandwidth ever observed 
and thus constitutes a bright source of quantum-correlated 
photons. 



II. EXPERIMENT 

At present the most successful squeezed light sources are 
based on parametric processes in nonlinear optical materials. 
In the regime of pulsed laser radiation the optical Kerr effect 
has recently enabled the observation of nearly 7 dB squeez- 
ing [24]. In the regime of continuous-wave laser radiation, 
squeezing of vacuum fluctuations by up to 10 dB was ob- 
served utihzing optical parametric oscillation (OPO) below 
threshold [25, 26]. 

The squeezed vacuum states investigated here were produced 
via type I OPO below threshold in a nonlinear standing wave 
cavity, similar to the experimental set-up in [26]. In order to 
observe high levels of squeezing the entire experiment was 
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optimized for low phase fluctuations and low optical loss. 
A schematic diagram of the experimental setup is shown in 
Fig. 1(a). The main light source was a 2 Watt continuous wave 
laser at 1064 nm. It was used to generate a second harmonic 
light field to pump the parametric squeezed light source and 
to provide the local oscillator (LO) beam for balanced homo- 
dyne detection. Both laser fields were sent through travelling 
wave filter cavities that were kept on resonance using Pound- 
Drever-Hall modulation-demodulation control schemes. The 
filter cavities therefore acted as optical low-pass filters that re- 
duced potential phase noise and provided clean gaussian fun- 
damental modes. The filtered second-harmonic pump beam 
had a power of up to 600 mW and was subsequently injected 
into the OPO cavity. 

In order to minimize the internal optical loss of the squeezed 
light source we designed the OPO as a monolithic (bi- 
convex) cavity, similar to the resonator used in [36]. It was 
made from a 7% doped Mg:0 LiNbOa crystal measuring 
6.5x2x2.5 mm"^ in dimensions with end faces of 8 mm radii 
of curvature (see inlay in Fig. 1). While the back of the crys- 
tal had a high reflection coating for both the wavelengths in- 
volved, we chose a rather low power reflectivity of 88 % at 
1064 nm (and < 1 % at 532 nm) for the coupling front face 
in order to increase the ratio of out-coupling rate and cavity 
round trip loss. The crystal was controlled to keep its phase 
matching temperature close to 60°C using peltier elements, 
and the main laser frequency was tuned to match the reso- 
nance condition of the OPO cavity. An assembly drawing of 
the crystal mounting is shown in Fig. 1(b). 
No control fields at or near the fundamental wavelength of 
1064 nm were injected into the squeezed light source during 
data taking. Thus all photons around 1064 nm leaving the 
OPO cavity originated from the parametric down-conversion 
process. The out-coupled squeezed states, i.e. the photons 
at 1064 nm, were separated from the 532 nm pump upon re- 
flection at a dichroic beam splitter (DBS) and subsequently 
analyzed by means of balanced homodyne detection. To this 
end the bright and spatially filtered LO beam was precisely ad- 
justed to the spatial mode of the squeezed beam. We achieved 
a fringe visibility of 99.8 % at the homodyne detector's 50/50 
beam splitter, which proves that the squeezed states were in- 
deed in a well-defined gaussian fundamental mode. 
A homodyne detector is the only known means to observe and 
characterize strongly squeezed states. Todays single photon 
detectors do not reach the same high quantum efficiencies and 
have an insufficient photon number resolution. It has been 
shown that homodyne detection data of quantum states can in- 
deed provide their /mW quantum information [31, 32, 33]. The 
phase angle 9 between the local oscillator (LO) and the signal 
beam defines an arbitrary quadrature of the signal's quantum 
state X{9) = cos{9)Xi+sin{6)X2 [21]. In our experiment 9 
was controlled with a piezo actuated mirror in the local oscil- 
lator path which allowed us to precisely set the phase angle. 
For example, to measure the squeezed quadrature X{Q) = Xi 
or the anti-squeezed quadrature X{7r/2) = X2- The collected 
data was readily converted into variances by a spectrum ana- 
lyzer that was also used to select a certain Fourier frequency 
/ and resolution bandwidth A/. The vacuum noise reference 
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FIG. 1: Experimental apparatus, (a). Schematic diagram of the ex- 
perimental setup. A small part of the 1064 nm output of the 2 W laser 
system was used as local oscillator while the majority was sent to an 
external second harmonic generation (SHG) cavity which provided 
the 532 nm pump beam to drive the OPO. Travelling wave filter cav- 
ities (FC1064 and FC532) in both paths were used to provide stable 
and well defined laser beams. Squeezed vacuum states at 1064 nm 
were generated by type I optical parametric oscillation (OPO) be- 
low threshold in a nonlinear monolithic cavity. A balanced homo- 
dyne detector was used to measure the states' quadrature variances. 
DBS: dichroic beam splitter; PD: photo diode; SA: spectrum ana- 
lyzer, (b), Exploded assembly drawing of the oven enclosing the 
squeezed light source. The nonlinear crystal, copper plates, peltier 
elements and thermal insulations are shown. The inlay shows a 
photograph of the monolithic squeezed light source made from 7 % 
doped MgOiLiNbOa. 



was measured with the same settings but with the squeezed 
field input blocked. 



A. Observation of high-purity strongly squeezed states 

The traces in Fig. 2(a)-(c) represent three different squeezed 
vacuum states from our OPO source in terms of pairs of 
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FIG. 2: Quadrature noise variances as measured by balanced ho- 
modyne detection. All traces were recorded at a Fourier frequency 
of /=5MHz, with a resolution bandwidth of A/ =100 kHz and 
video bandwidth of 100 Hz, and were normalized to the vacuum 
noise. The data still include electronic detector dark noise, hence all 
traces represent direct observations without any data post-processing, 
(a), Squeezed vacuum state that was generated with a second har- 
monic pump power of Psh=30 mW. A nonclassical noise reduction of 
2.9 dB below vacuum noise was observed, with the respective anti- 
squeezing of 2.9 dB above shot noise, (b), Squeezed vacuum state 
generated with an increased pump power of Psh=150 mW. Squeezing 
of 6.2 dB was measured with 6.7 dB of anti-squeezing. (c). Driving 
the OPO at maximum pump power Psh=600mW yielded squeezing 
and anti-squeezing of 1 1 .5 dB and 1 6.0 dB respectively. The subtrac- 
tion of the detector dark noise would further increase the value of 
maximum squeezing to 11.6dB. 



squeezed (bottom trace) and anti-squeezed (top trace) vac- 
uum normalized quadrature variances Vi = A^Xi/A^Xi vac 
and V2 — A^X2/A^X2,vac, respectively. For a pump 
power of 30 mW (Fig. 2(a)) the squeezed variance (Vi) 
was 2.9(±0.1)dB below and the anti-squeezed variance (V2) 
2.9(±0. l)dB above the vacuum noise level. Within the mea- 
surement error bars of ±0. 1 dB the state shows minimum 
uncertainty and therefore perfect purity. Upon driving the 
OPO with 150 mW pump power (Fig. 2(b)) we observed 
normalized variances of -6.2(±0.1)dB and 6.7(±0. l)dB, re- 
spectively. The strongest and unprecedented squeezing level 
was measured with a pump power of approximately 600 mW 
(Fig. 2(c)). With this setting 16.0(±0. l)dB of anti-squeezing 
and a noise reduction of 11.5(±0. l)dB below the vacuum 
noise level was achieved. 

In practice, squeezed vacuum states generally cannot have 
perfect purity due to both residual optical loss and residual 
phase noise. The former results in an admixture of the vac- 
uum state and the latter results in a mixture of phase space 
rotated states [25, 27]. Here, we make the attempt to describe 
all the three squeezed states measured by a fixed amount of 



optical loss assuming zero phase noise. This approach is rea- 
sonable for the following reasons, i) previous investigations 
have shown that our set-up is stable and phase noise is low 
[26], ii) in contrast to squeezing in the audio-band [29], at 
MHz Fourier frequencies no displaced or thermal states are 
present that could mix into the state, and iii) no evidence was 
found for a non-Gaussian statistic in our measured data. All 
3 pairs of variances observed can then be consistently calcu- 
lated using Vi = rjjV^ + (1 — 777), where V( = I/V2 repre- 
sent the pure state variances before any decoherence has oc- 
cuiTed. Here, (1 — 777) represents the overall contribution of 
the vacuum state, in power. The detection efficiency 77 ac- 
counts for all losses outside the squeezed light source and 7 
is the OPO cavity escape efficiency accounting for all losses 
inside the squeezed light source. We find that all of the mea- 
sured variances (shown in Fig. 2) can indeed be modelled by a 
vacuum noise contribution as low as (1—777) = (4.8 ±0.2) %, 
independent of the pump power used. This result supports 
our assumption that phase noise is not significant in our set- 
up. The loss value of 4.8% is consistent with an indepen- 
dent loss analysis of all the optical components used in our 
set-up, however, an independent measurement of the quantum 
efficiency of the photo diodes of the homodyne detector typ- 
ically has large uncertainties, and we estimate this value to 
(98 ± 2)%. Such a low optical loss should enable the obser- 
vation of more than 13 dB of squeezing, however, pump pow- 
ers above 600 mW stably controlled and within a clean funda- 
mental mode were not available in our experiment. Figure 3 
presents the Wigner function [30, 31, 32] of the squeezed vac- 
uum state with variances of -11.5dB and 16dB with respect 
to the vacuum noise variance. The Wigner function represents 




FIG. 3: Wigner function of the squeezed vacuum state produced 
by our OPO. The projections (filled curves) onto the two quadratures 
yield the gaussian probability distributions with variances of - 1 1 .5 dB 
and 16 dB relative to the projections belonging to a pure vacuum state 
(dotted curves). 
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one of several possible phase-space probability distributions 
for quantum states. Since it is Gaussian for vacuum-squeezed 
states, it is particularly easy to obtain simply by measuring its 
second moments. The entire statistical properties of the state 
can than be summarized in a finite dimensional space, namely, 
by the covariance matrix. Hence, for squeezed vacuum states 
mixed with the ordinary vacuum state only the variances of the 
squeezed and anti-squeezed quadrature are required to calcu- 
late the Wigner function which is then given by 



: exp 
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where A'^Xi are again the variances of the quadratures Xi 
and X2, respectively, with the vacuum noise variance of 1/4. 
The vertical axis in Fig. 3 represents the quasi-probability den- 
sity of the quadrature phase space. The Wigner function is a 
phase-space ^MaiZ-probability distribution, since the quadra- 
tures cannot simultaneously have precisely defined values. 
What can be precisely measured is one projection of the 
Wigner function onto the axis of the phase space at a time. 
The projections are the probability distributions of the particu- 
lar quadrature, and here, have a Gaussian shape with squeezed 
or anti-squeezed variances with respect to the vacuum noise 
variance (dotted curves). Contour lines of the Wigner function 
have an elliptical shape and are thus squeezed in comparison 
with the circular contour lines of the vacuum state (not shown 
in Fig. 3). 



B. Squeezed vacuum states in the photon number basis 

All photons that are present in a squeezed vacuum state pro- 
duced by OPO stem from parametric down conversion of sec- 
ond harmonic pump photons. This is also true for a mixture of 
a squeezed vacuum state with the ordinary vacuum state. For 
strongly squeezed states not just the detection of photon pairs 
but also of higher-order photon pairs is expected. However, 
the overall loss of photons must be sufficiently low to enable 
such an observation. Since the Wigner function represents the 
full quantum information of a state it can be used to calculate 
the state's density matrix p„_,„ in the photon number (Fock) 
basis, which is another complete description of the state. The 
photon number distribution P(n) can be obtained from the 
diagonal elements of the density matrix. Figure 4 shows a 
plot of the first matrix elements (absolute values) up to pho- 
ton numbers n,m ^ 10 of the observed squeezed state with 
noise variances of -1 1.5 dB and 16 dB relative to the vacuum 
noise. In order to calculate the density matrix [21, 34] photon 
numbers of up to n = m = 170 were considered, which en- 
sured convergence of all the matrix elements. Similarly, the 
density matrices for the additional two states have also been 
generated. Their diagonal elements, i.e. the photon number 
distributions of the states, are depicted in Fig. 5. The for- 
mula used and the numerical values of the density matrix el- 
ements are given in the supplementary information. For the 
1 1 .5 dB squeezed state (red) a strong odd/even photon num- 
ber oscillation is visible. Even the probability of detecting 




FIG. 4: Density matrix of the squeezed state of Fig. 2(c). Plot of the 
absolute values of the density matrix elements with photon numbers 
up to 10 for the -11.5 dB squeezed and 16dB anti-squeezed state. 
The red diagonal elements give the probabilities of finding 0,1,2. .. n 
photons. Due to the state's high purity, the probabilities for finding 
an odd number of photons is significantly suppressed. 



n = 10 photons is significantly higher than the probability 
of detecting n = 9 photons. This oscillation further contin- 
ues for greater photon numbers but is not shown here. For 
our squeezed states with less squeezing the oscillation flat- 
tens out at lower photon numbers, but the odd/even oscillation 
is even more pronounced at lower photon numbers. For the 
2.9 dB squeezed state the probability of finding n — 1 photon 
is as small as 0.5%. Note that for displaced squeezed states 
additional photons are added to the state and the odd photon 
numbers can not all have zero probability, even if these state 
were pure. However, weakly displaced squeezed states show 
photon number oscillations with larger periods. Such states 
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FIG. 5: Photon number oscihation. The photon number distribu- 
tion of the three squeezed vacuum states that were generated in our 
setup clearly show oscillatory behavior. With increasing squeezing 
the number of higher order photon pairs grows. For clarity the inset 
shows a zoomed view restricted to probabilities up to 0.1. 
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were originally theoretically analyzed in [35]. 
The photon numbers calculated here are detected per band- 
width in Hertz and per detection time interval in seconds. Ob- 
viously, the higher the squeezing bandwidth the higher the 
number of photons per second. 



in [28, 29]. 

The spectrum of the squeezed and anti-squeezed quadrature 
variances (Vi and V2, respectively) for the below-threshold 
OPO can be calculated by [25] 



C. Observation of high-bandwidth squeezing 

For applications in high speed quantum communication and 
information processing high photon counting rates or, equiva- 
lently, fields with squeezing over a broad Fourier spectrum are 
essential. The half width half maximum (HWHM) squeezing 
bandwidth can be defined as the spectral width over which the 
squeezed variance increases from its lowest variance at small 
Fourier frequencies Vi to Vi+0.5-(l — Vi). Here the variances 
are again normalized to unity vacuum noise. Although contin- 
uous variable quantum states are often quoted to provide the 
possibility of high bandwidths, only modest values of the or- 
der of a few tens of MHz have been realized, e.g. in [36, 37]. 
Figure 6 shows the spectral analysis of the strongly squeezed 
field from our source with a HWHM squeezing bandwidth of 
170 MHz. In contrast to Fig. 2 here the electronic dark noise 
of our homodyne detector was subtracted. At ^100 MHz the 
dark noise reached the measured quantum noise. At frequen- 
cies below 5 MHz the measurement was limited by the reso- 
lution bandwidth chosen. Since our set-up did not use control 
beams which might introduce technical noise at low frequen- 
cies, a white squeezing spectrum can be expected for lower 
frequencies if parasitic interferences with frequency shifted 
scattered photons are sufficiently suppressed as demonstrated 
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FIG. 6: High bandwidth squeezing spectrum. Squeezing (bottom 
trace) and anti-squeezing (top trace) are shown relative to the vac- 
uum noise variance. The measurements were performed from 5 MHz 
to 100 MHz. This was limited by the fact that the dark noise clear- 
ance of the homodyne detector was too low towards higher frequen- 
cies. The traces are averages of 3 measurements each done with a 
resolution bandwidth of 1 MHz and video bandwidth of 3 kHz. The 
data was fitted to a model and revealed a squeezing bandwidth of 
170 MHz (thin line). 
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where P532 is the power of the second harmonic pump and 
Pth is the pump power needed to reach the OPO threshold. 
The detection efficiency 77 is calculated as the product of the 
propagation efficiency, the square of the homodyne visibility, 
and the quantum efficiency of the homodyne detector photo 
diodes. The OPO escape efficiency 7 = T/{T + L) is de- 
termined by the transmittance T of the out coupling mirror 
and the intra-cavity round trip loss L. The dimensionless pa- 
rameter K ~ 27r/ /k is the ratio of the Fourier frequency / 
and the decay rate k, = {T + L)c/l of the cavity with round 
trip length I. The theoretical curves in Fig. 6 were modeled 
using equation (2), with the best agreement obtained upon as- 
sumption of the following values: a pump beam at 53.5% of 
the threshold power, the output coupler transmittance of 12%, 
the round-trip loss within the crystal of approximately 0.1%, 
and a total detection efficiency r/7=95.2% in agreement with 
our analysis presented in section II A. The squeezing band- 
width of our source is given by the Fourier frequency where 
the squeezed variance has increased to -2.7 dB. According to 
the theoretical curve in Fig. 6 this bandwidth is 170 MHz. This 
result can be used to calculate the rate of down-converted pho- 
tons from a full longitudinal mode of our source. To this 
end, we split the half free spectral range of our source of 
5.5 GHz into 55000 frequency bins of 100 kHz each. For ev- 
ery bin z, the photon number distribution Pi{n) was calcu- 
lated. By summing these distributions, we find that the spec- 
trally weighted representative photon number distribution of 
the central OPO cavity mode around half the pump light fre- 
quency Pf{n) = 1/55000 • Pi{n). This distribution can 
now be used to calculate the mean number of photons per sec- 
ond and bandwidth averaged over the band of 5.5 GHz. Multi- 
plying with 5.5 GHz provides the rate of down-converted pho- 
tons from this mode of as high as 2.79 • lO^s^^, corresponding 
to a power of 52 pW. 

Another relevant value is given by the mean number of cor- 
related photons in the case of a single photon detection. For 
a parametric down-conversion source operated far below its 
threshold and with zero photon loss this number is two be- 
cause the detection of a single photon proves the presence of, 
in total, two photons. In order to calculate this number for 
our source we set P/(0) — and calculate the mean value 
of the renormalized remaining distribution. For the 1 1 .5 dB 
squeezed state presented in Fig. 6 we obtain (n) / = 5.93. Im- 
plying that whenever a single photon detector clicks (possibly 
without any photon number resolution) this event proves the 
presence of, on average, almost 6 quantum correlated photons, 
including the detected photon(s). 
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in. DISCUSSION 

In order to observe strong squeezing levels, we operated an 
optical parametric oscillator comparatively close to its thresh- 
old. In this regime not just single photon pairs, but a signif- 
icant number of multiple or higher order photon pairs were 
produced. This contrasts typical parametric down-conversion 
experiments in the photon counting regime where the prob- 
ability of observing more than two quantum correlated pho- 
tons is negligible. The density matrix in the photon number 
basis of our - 1 1 .5 dB squeezed state not only reveals strong 
photon number oscillations but also quantum correlations be- 
tween photons in the following way: Whenever a photon was 
photo-electrically detected by the homodyne detector, there 
was an almost unity probability that an additional, odd number 
of photons were physically realized and were detected by the 
same detector In photon counting experiments our squeezed 
state might therefore be used to overcome the problem of 
insufficient photon counter efficiencies. The unprecedented 
quality of our squeezed states will also enable a great step for- 
ward in many applications currently under discussion. In high 
power laser interferometers for gravitational wave detection 
[4, 5, 6] the application of an -1 1.5 dB squeezed state reduces 
the shot noise equivalent to a laser power increase by a factor 
14. The interference of two such states on a 50/50 beamsplit- 
ter would produce strong Einstein-Podolsky-Rosen entangle- 
ment [12, 13, 38] and its high bandwidth may result in the 
ultra-fast generation of quantum keys [39]. Our squeezed light 
source may also find applications in future hybrid quantum 
information networks. Its well-defined fundamental Gaussian 
spatial mode yielded a high fringe contrast with a spatially fil- 
tered local oscillator of 99.8% and enables efficient coupling 
to other systems, for example to atoms coupled to optical cavi- 
ties [40, 41, 42], in order to realize quantum interfaces of light 
and matter. 
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APPENDIX 

In order to calculate the density matrix elements Pm.n we 
followed the prescription given in [34], whereby the entries 
can be obtained using the following equation: 



jmW. \ - (_y)[(™-n)/2]+2a[/n-2Q 

^"'""YVV^^^ a!(„_2a)!(a+=^)! ' 
where 14 = + 1/2, Vp ^ Vp + 1/2, [/ = 1 - - 

and T = Note, that Vx and Vp correspond to the 

linear quadrature variances of the quantum state as introduced 
in equation (2), but with the vacuum variance set to 1/2. This 
expression for the density matrix in the Fock state basis can 
be derived from the covariance matrix by first noting that the 
Husimi Q-function in the coherent state basis is the generating 
function for the density matrix in the Fock state basis. The co- 
variance matrix needs to be transformed from the quadrature 
basis into the coherent basis to obtain the above formula. 
We calculated pm.n for the three states in Fig.2a-c. For the 
first state (Fig. 2a) we chose vacuum normalized variances of - 
2.84 dB and 2.94 dB being within our experimental error bars. 
These numbers correspond to an originally pure state with 
3.05 dB squeezing and 3.05 dB anti-squeezing that was de- 
graded by 4.8% of optical loss. For the remaining two states 
we chose the variances according to the values in Fig. 2b (- 
6.2dB and 6.7dB) and Fig.2c (-11.5 dB and 16.0dB). The 
calculation was done with photon numbers up to 170 to en- 
sure the convergence of all entries, but the matrices presented 
below were truncated at n=m=10. For clarity the diagonal ele- 
ments representing the photon number distribution (as shown 
in Fig. 5) are marked in red. 
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■ .7538 -.3360 .1834 -.1056 .0622 -.0372" 

.0131 -.0101 .0071 -.0048 .0032 

-.3360 .1500 -.0820 .0473 -.0279 .0167 

-.0101 .0078 -.0055 .0037 -.0025 

.1834 -.0820 .0449 -.0259 .0153 -.0092 

.0071 -.0055 .0039 -.0026 .0018 

-.1056 .0473 -.0259 .0150 -.0089 .0053 

-.0048 .0037 -.0026 .0018 -.0012 

.0622 -.0279 .0153 -.0089 .0053 -.0032 

.0032 -.0025 .0018 -.0012 .0008 

-.0372 .0167 -.0092 .0053 -.0032 .0019 _ 

" .3026 -.1946 .1532 -.1272 .1082 -.0933' 

.0126 -.0140 .0143 -.0140 .0135 

-.1946 .1256 -.0993 .0828 -.0707 .0613 

-.0140 .0156 -.0159 .0157 -.0151 

.1532 -.0993 .0789 -.0660 .0566 -.0493 

.0143 -.0159 .0162 -.0160 .0155 

-.1272 .0828 -.0660 .0555 -.0478 .0417 

-.0140 .0157 -.0160 .0158 -.0153 

.1082 -.0707 .0566 -.0478 .0413 -.0362 

.0135 -.0151 .0155 -.0153 .0148 

-.0933 .0613 -.0493 .0417 -.0362 .0318 
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